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Abstract. The low temperature properties of fermionic solids are governed by 
spin exchanges. Near the melting transition, the spin-exchange energy increases 
as well as the relative contribution of large loops. In this paper, we check the 
convergence of the multi-spin exchange model and the validity of the Thouless 
theory near the melting of the Wigner crystal in two dimensions. Exchange 
energies are computed using Path Integral Monte Carlo for loop sizes up to 8, 
at r s = 40, 50 and 75. The data are then fitted to a geometric model. Then 
the exchange energies are extrapolated to larger sizes in order to evaluate their 
contributions to the leading term of the magnetic susceptibility and the specific 
heat. These results are used to check the convergence of the multi-spin exchange 
model. 



1. Introduction 

3 He atoms as well as electrons carry a spin 1/2. In their solid phases, exchange 
processes are responsible for the thermodynamics at low temperature. The Thouless 
theory explains how a Multi Spin Exchange (MSE) effective Hamiltonian results from 
these spin exchanges pQ. The applicability of such an approach is based on the very 
different scales of the exchange energies and the phonon energies (at least 3 orders 
of magnitude). In the strong coupling limit (r s — > oo for the Wigner crystal), a 
semiclassical calculation (WKB) predicts such very small exchange energies. Moreover, 
in this limit, only the two or three body exchanges dominate depending on the 
dimension and the lattice type, leading to the standard Heisenberg model (for spin 1/2, 
three body exchanges are equivalent to pair exchanges). Then, at weaker coupling, 
other loop exchanges arise that may lead to a magnetic phase transition. In general, 
more and more exchanges are important near the melting transition. The problem 
we address here is how the MSE Hamiltonian evolves as we approach the melting 
transition and whether there is a precursor to melting in the magnetic Hamiltonian. 

Exchange energies have been obtained for the bulk 3D 3 He[3] and for 3 He 
adsorbed on graphite^. The exchange energies have been obtained as well for the 2D 
electron Wigner crystal|5]-[H|- F° r most densities, from the semiclassical limit to near 
melting, the WKB and Path Integral Monte Carlo (PIMC) evaluations assume a stable 
solid of particles (a triangular lattice in 2D) without taking into account the Fermi 
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statistics. For these densities, only small exchange loops with 6 or fewer particles 
have significant energies. These calculations determine the exponential dependence of 
exchange energies versus the density; the longer loops decreasing more rapidly in the 
strong coupling limit. The MSE model explains the main features of the magnetic 
transition of the second solid layer of helium 3 adsorbed on graphite[7j. 

One problem that has been encountered near melting is that the PIMC 
calculations become increasingly difficult to converge. This is because the 
Boltzmannon solid melts before the fermion solid. Adding Fermi statistics stabilizes 
the solid in a wider density region. (Or alternatively, imposing antisymmetry increases 
the energy of the liquid phase into which the crystal melts.) In order to stabilize 
the crystal, we assume a ferromagnetic spin polarization of "spectator" atoms (non- 
exchanging atoms) and enforce antisymmetry on those paths. When comparison is 
possible, only small changes of exchange energies are found with this procedure. The 
largest differences are seen at r s = 50 for the system without statistics where the solid 
is close to melting (melting of the bosonic wigner crystal occurs at r s w 60). The 
restriction on the path resulting from Fermi statistics does not modify the exchange 
energies very much. The statistics do not affect the local correlated motions of particles 
in the vicinity of the exchanging ones. But statistics have long range effects that 
stabilize the solid and forestall the melting. 

The next question is whether the Thouless theory can be applied near the melting 
transition. Indeed, a possible scenario could be that the energy difference associated 
with a given spin-exchange, say Jp, becomes comparable to the zero-point motion 
kinetic energy (ki ' 1 ). In that complicated coupling between spin states and 

phonon excitations could arise. Instead our numerical results demonstrate, that all 
the exchange energies are much smaller than ki°^ even at melting. Even if the 
exchange energy increases exponentially with the inverse density (WKB predicts 
— log(J) cx r^ 2 ), each value stays well below ki°\ Thus, up to the melting transition, 
each exchange can be seen as an independent event. 

Now, near the melting transition, the number of important exchanges increases 
considerably. In this paper, we address the question of the convergence of this MSE 
model versus the number of exchanges. Indeed, if the exchange energies decrease 
roughly exponentially with the loop size n, the number of loops of size n does increase 
exponentially. The convergence of the MSE model is thus questionable. Does the 
total spin-exchange energy diverge when the crystal melts? Is the existence of long- 
exchanges a mechanism for the melting transition or a by-product? We try to answer 
these questions from the analysis of the two-dimensional Wigner crystal exchange data. 
In section II, we present the results of Path Integral Monte Carlo (PIMC) simulations 
and we give the formula to test the convergence of the MSE model. In section III, 
different models are used to fit PIMC results in order to extrapolate to larger loop 
sizes. 

2. PIMC Results and Thermodynamics 

For a given permutation P, the exchange energy Jp calculated from PIMC needs to be 
monitored carefully with respect to several parameters such as the number of particles 
and the time step r. There are additional parameters associated with the Bennett 
method. Fortunately, those parameters change smoothly with the density or the type 
of exchange or even the system (Helium atoms or electrons). The most important 
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effect is the (3 = 1/T dependence of the results. Indeed, the calculations are done 
at a temperature sufficiently low so that the phonon excitations are non-existent but 
sufficiently high that the spins have not magnetically ordered. In these conditions, the 
exchange energy should be independent of the temperature. However, near melting 
(the precise density depends on the statistics) results show a linear increase versus (3 
(note that each independent run at a given temperature is well converged). This is in 
fact interpreted as a signature of melting that would most probably occur for larger 
systems. The following tables contain only exchange values converged and independent 
of p. 

Table □ shows the PIMC results. Shown are the energies versus exchange for the 
3 densities considered. It is seen that exchange energies decrease as r s increase or at 
fixed r s , when the loop size increases. The strong dependence on the shape of the 
exchange we examine next. 

The overall importance of the exchange processes can be measured by the leading 
term of the high temperature (HT) expansion of the magnetic susceptibility (the 
Curie- Weiss temperature 9) and the specific heat: 

where 

Jn ' = ^ t n shapeJn P ^ (2) 

shape of n sites 

with n s hape is the number of loops per site of a given shape and j^ sha P e "> ^ ne associated 
exchange value. 

The leading HT-term of the specific heat is defined as 

^gw„*(JaA a (3) 

t^oo Nk B 4 V T J w 
where Jc v is given by a positive quadratic form of the exchange values: 

■l<\ }2 J - }l -r J r (4) 

Unfortunately, no simple formula exists for the matrix M. At very large r s (WKB 
limit) only Ji and mainly J3 dominate in a ferromagnetic phase and \9\ — Jc v = 
J2 — 2J3. At r s ~ 200, J4 is large enough and the system becomes anti ferromagnetic. 



3. Fits of exchange energies to analytic formula 

In order to determine the contribution of larger exchanges to the magnetic 
susceptibility and specific heat in Ecis lll3l we want to understand how exchange 
energies vary with the loop size and shape. 

Exchange energies depend mainly on the loop size n. Then, for a given loop size 
n, the exchange is much larger if the exchange contains only smooth angles, that is 
when its area A is maximum (^4 = (n — 2 + 2p)Ao, where Ao is the area of the enclosed 
triangles and p is the number of lattice sites inside the loop). In order to quantify this 
shape dependence we assume that each local section of the exchange contributes to 
the exchange probability and use the following formula: 

log(J^*) = a + a n n + a p p s + ^ a vN s ,v (5) 



Multi-spin exchange model the near melting transition of the 2D Wigner crystal 4 



J 


mult 


r s = 


40 


r s = 


50 


r s = 


75 


2 


3 


1402. 


2% 


323. 


3% 


13.4 


3% 


3 


4 


1184. 


3% 


281. 


3% 


14.2 


2% 


4 


6 


809. 


3% 


165. 


5% 


6.2 


3% 


5 


12 


295. 


6% 


42. 


7% 


1.1 


6% 


6i 


2 


476. 


5% 


88. 


9% 


2.3 


11% 


6 2 


12 


143. 


7% 


17.1 


8% 


0.24 


8% 


63 


12 


92. 


8% 


9. 


9% 


0.116 


9% 


64 


4 


40. 


13% 


3.2 


8% 


0.021 


8% 


7i 


12 


105. 


9% 


14.4 


6% 


0.16 


18% 


7 2 


12 


72. 


11% 


7.8 


9% 


0.076 


21% 


7 3 


24 


43. 


11% 


3.6 


9% 


0.018 


13% 


7 4 


12 


37. 


21% 


2.8 


16% 


0.0037 


19% 


7 5 


24 


20. 


15% 


0.97 


11% 


0.004 


17% 


81 


6 


130. 


11% 


15.3 


7% 






82 


6 


29. 


13% 


1.74 


9% 






83 


12 


19. 


12% 


1.19 


10% 






84 


12 


23. 


10% 


1.72 


7% 






85 


24 


40. 


13% 


3.1 


9% 






8 6 


12 


44. 


12% 


3.5 


11% 






87 


6 


1.3 


21% 


0.058 


17% 






8 8 


12 


4.3 


26% 


0.32 


26% 






89 


24 


6.4 


24% 


0.25 


21% 






810 


12 


23. 


11% 


1.5 


14% 






811 


24 


7.6 


14% 


0.25 


20% 






812 


24 


18. 


18% 


0.97 


11% 






813 


12 


8.8 


16% 


0.24 


15% 






814 


24 


22. 


10% 


2.5 


8% 






815 


24 


9.8 


14% 


0.51 


10% 






816 


12 


14. 


12% 


1.08 


8% 






9i 


4 


75. 


15% 


5.2 


10% 





Table 1. Exchanges energies (in 10 9 Ry) obtained by PIMC simulations for the 
fully polarized 2D Wigner crystal. Results are for 64 electrons, and respectively 
P = 1/T = 10000, 16600, 20000 Ry- 1 , time step r =40, 66, 100 Ry- 1 for r s = 40, 
50, 75. The multiplicity is the number permutation P of the same shape per site, 
including p- 1 symmetry. All exchanges are for nearest neighbor on the triangular 
lattice. For 2 through 6 particle exchange there is only one sort of diagram, but for 
larger exchanges there are multiple possibilities which we have labelled arbitrarily 
with an index. The authors can be contacted for further details. 



Here p s is the number of non-exchanging sites contained inside the loop of shape s 
and N s _ v is the number of patterns of type v encountered in the loop of shape s. The 
a's are parameters determined from a least-squares fit minimizing: 




where Nj are the number of degrees of freedom (the number of exchanges minus the 
number of fitting parameters), and a n , s is the Monte Carlo error of log(J„ iS ). We now 
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Table 2. Top Table : Fit of exchange energies according to Eqs. 15161 a n is the 
main varying parameter. Uncertainties are in parenthesis. Middle Table : Same 
for fits with the 2 links patterns. For each r s , the third column is the ratio Vi/vo- 
within error bars, they are independent of r s for this model. Bottom Table : same 
for 3 links patterns. For each r s , the third column is the ratio Wi/wo, the ratio 
depends on r s for patterns with smooth angles, and not for patterns with sharp 
angles. 



discuss the results for the fit allowing for 3 successively more detailed patterns. 

Fitting without patterns. Without using any of the parameters that depend 
on the shape of the exchange (a v = 0), only 3 parameters remain: ao, a n and a p . 
The results are given in table [3 Top. Such a fit is not accurate as can be seen by the 
large value of x 2 - As an interpretation of the parameters, exp(ao) is the prefactor 
for the initiation of an exchange, a n is the action needed for adding a link onto the 
exchange, and a p , is a term which favors non-compact exchanges. 

Fitting with 2-links patterns. The simplest patterns v consist of the angles 
subtended by 2 adjacent links. On a triangular lattice there are three possible cases, 
since two adjacent links can have angles of 0, 60, 120 or 180 degrees (see FigQ but 
the parameter v±so only occurs in Ji . Also a n is becomes redundant because it can be 
determined from the angles: N StV = n, the number of angles equals the number 
of sites in the loop. Thus we are left with five parameters. The results for this fit are 
shown in table [^-Middle. 

The mean error is significantly reduced by including the angles in the fit. We 
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Figure 1. List of the three possible patterns made of two links (last line) and 
the eight patterns made of three links. 



find that «o depends weakly on density, a p is not very significant. Note that the ratio 
Vi I vq is almost independent of r s . 

Fitting 3-links patterns. In our most elaborate fit, the possible configurations 
of three adjacent link are used (see Fig. The triangular pattern w 8 occurs only in 
the 3-body loop, thus is removed from the fitting procedure. Note that if the patterns 
for 3-links are used, those of 2-links are redundant. Here again a n is a redundant 
parameter. Hence, we are left with 10 parameters for 28 exchanges for which we have 
computed exchange frequencies. At r s — 75, we have not computed enough different 
exchanges, to do this fit so we only discuss it at smaller values of r s . 

The results are shown in Table |3-Bottom. This fit is now reasonable as the mean 
standard deviation is of order of one, i.e. within statistical error bars of the PIMC 
results. Clearly the parameter a p can be removed. In this fit, the ratios Wi/wo depend 
on r s except for W3, 105, w-j. 

4. Expanding to larger loop cycles 

From the fits of the previous section, it is now possible to evaluate, at high 
temperature, the leading term 9 of the magnetic susceptibility (see EqQJ and Jc v 
of the specific heat (see Eq|3J) by including much larger loops. Those terms give an 
estimate of the temperature of the spins begin to become correlated (note that the 
expected ground state is a spin liquid) . We test here if the MSE model still converges 
when more and more exchanges are involved as it appears near the melting transition. 

We have calculated all self avoiding walks of length up to size n = 22 on the 
triangular lattice, with their shape properties (p, angles v and w). 

Fig[5}left shows the sum jj of all energies of loops of size n. jj decreases at 
large r s , but clearly increases exponentially at r s = 40. Nevertheless the Curie- Weiss 
temperature seems to converge at all densities (Fig|2}Middle) even if strong oscillations 
appear at the r s =40. 

The leading term of the specific heat Jc v (Eq|2J can be used to scale 
energies/temperatures. A divergence of this term means a huge energy associated with 
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Figure 2. Left : Total Exchange energy of loops of length n versus n (see Eg 151. 
Middle : Contribution of loops of length n to the Curie- Weiss temperature versus 
n (see Eq0. Right: Contribution of loops of length n to Jq v versus n (see Easl3t 
11*1 . The full lines stand for the model with 3-link patterns, dashed lines for the 
model with 2-link patterns, and the dotted line for the model without patterns 
(see text). 



exchanges, and likely the melting of the solid. Fig^right shows a good convergence of 
Jc v again at large r s . The 2-links and 3-links models provide rather different results. 
This means that Jq v is a rather subtle combination of exchange energies. Still all 
models seem to converge at r s — 50. At the smallest r s — 40, the 2-links model 
diverges, whereas the 3-links model seems to converge though with strong oscillations. 

5. Conclusion 

Exchange energies can be computed with Path Integral Monte Carlo up to the melting 
transition. Use of Fermi statistics is necessary to stabilize the solid (non-local effect), 
but it does not change the exchange energies (local contribution). As one approaches 
the melting transition, an increasing number of different exchanges becomes impor- 
tant. The Thouless theory allows one to write a multi-spin exchange model to describe 
the low energy physics. But the convergence of this model is not guaranteed and must 
be verified. We have checked that for the 2D Wigner crystal, the Curie- Weiss tem- 
perature 9 and the first-order contribution to the specific heat (Trace(i/— < H >) 2 ) 
converge for the lowest densities. The large oscillations seen both these quantities at 
r s = 40 hint at a possible divergence in the thermodynamic properties at higher den- 
sities. Though the MSE model converges at lower densities, a divergence of the model 
is possible at the critical melting density, consistent with the picture that multi-spin 
exchange is itself a mechanism for melting. 
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